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ADAPTIVE  TRACKING  ALGORITHM  FOR  TRACKING 
AIR  TARGETS  WITH  SEARCH  RADARS 


INTRODUCTION 

Since  the  time  of  the  Wiener  Filtering  Theory,  which  was  based  on  stationary 
processes,  there  ha-i  been  increased  interest  in  adaptive  filters.  Some  of  the  better  known 
adaptive  systems  are  adaptive  antennas  [1],  phase-lock  loops  [2],  and  the  Kalman  filter 
[3] .  Recently  there  has  been  considerable  interest  in  adaptive  tracking  systems  [4-12] . 
Most  of  these  adaptive  systems  use  a  feedback  loop  to  adjust  the  parameters  in  the  sys¬ 
tem.  A  popular  method  of  adjusting  the  systems’  gains  is  to  use  a  least  mean  square 
(LMS)  error  criterion  and  minimize  it  with  respect  to  the  gains  [1,4] .  This  report  de¬ 
scribes  an  adaptive  a  -  0  filter  based  on  this  LMS  error  criterion.  Although  its  operation  is 
similar  in  principle  to  the  adaptive  antennas  [1]  and  tracking  system  [4] ,  there  are  marked 
differences.  We  begin  by  reviewing  the  a  -  (3  filter. 


REVIEW  OF  THE  a  -  0  FILTER 

The  a-B  filter  is  defined  as  [13,14] 


=  smoothed  position 
vt(k )  =  smoothed  velocity 

xp(k)  =  predicted  position 
xm(k)  =  measured  position 
*m(*)  =u(k)  +  u>(k) 


CD 


(2) 


Note:  Manuscript  tubmitted  July  26,  1974. 
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u(k)  -  true  position  (band-limited  but  unknown) 
w(k)  ~  zero  m!»an  white  Gaussian  measurement  noise 
T(k)  =  time  between  samples 
a,P  =  system  gains. 

Applying  the  z -transform  to  Eqs.  (1)  and  (2),  we  find  that  the  transfer  functions  of  the 
system  are 


—  _  —  _ _ 1= - - - 

x  xm(t)  ;2-Z(2-a-0)  +  (!-«) 


H  m  Mf) _ T(k)  Z{Z  " 1}  ( 

u  z2 -z(2-er-0)  +  (1 -a) 

The  transfer  functions,  Eqs.  (3)  and  (4),  are  placed  into  standard  notation  for  a  second 
order  system: 


H  =  - - - .  ({ 

{* }  z2  -  2ze_^wo'nt*)  cos  udT(k)  +  e~^onk) 

Equating  terms  in  the  denominators  of  Eqs.  (3)  and  (4)  with  that  of  Eq.  (5),  we  obtain 

a=l.e-^«oW  (( 

0  =  1  +  e-%<»onk)  „  2e^uoT(k)  cos  ^[k),  (' 

or  conversely, 


2 

/ ^  1 

+ 

s 

i 

H 

> 

1 _ » 

1 

rw-1  i 

2  -  a  -| 

2 
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where  £,  coa-,  and  a>0  are  the  classic  damping  coefficient,  damped  natural  frequency,  aid 
natural  frequency  of  a  second  order  system.  In  Ref.  13,  it  was  suggested  that  if  one  set 
a,  then  0  could  be  found  by 


0  = 


a2 


(2 -a) 


(ID 


Substituting  this  expression  into  Eq.  (8),  we  find  that  £  varies  from  0.707  to  0.86  as  a 
varies  from  zero  to  unity.  We  find  that  by  using  Eq.  (11),  the  system  remains  near  critical 
damping  for  all  values  of  a  and  that  a  and  T(k)  controls  the  system  bandwidth.  This 
concludes  the  review  of  the  a  -  0  filter. 


LEAST-MEAN-SQUARE  ERROR  CRITERION 


The  least  mean  square  error  criterion  which  is  defined  as  the  expected  value  of  the 
square  of  the  difference  between  the  target’s  predicted  and  measured  position  was  chosen 
as  the  system’s  preformance  measure; 


g=£{[xp(fc  +  l)-xm(fe  +  l))2}. 


(12) 


The  gradient  of  &  with  respect  to  the  system’s  gains  is  found  by  interchanging  the  ex¬ 
pected  value  and  derivative  operation; 


+  l)-xm(k  +  l)] 


3xJk  +  1)\ 

- -  =0 

3a 

(13) 

a*p(fe  +  in 

— c - •=  n 

30 

(14) 

The  gradients  are  set  to  zero  and  we  solve  for  the  system’s  gains  which  yield  minimum 
error.  If  Eq.  (1)  is  substituted  into  Eq.  (2),  xp(k  +  1)  becomes 


xp{k  *  1)  =  xp(fe)  +  T(k  +  l)vs(k  -1)4  «[xw  (fe)  -  xp(ft)] 
0T(k  +  l)[zm{k)  -  xPik)] 


T(k) 


(15) 


where  x  (k)  =  xt{k  -  1)  +  T(k')vf(k  -  1). 


The  partial  derivatives  of  xp(k  +  1)  are  proportional  to 


dx  (k  + 1)  dx  (k  1) 

■  — 1 . .  —  ■  II  — 


dct 


30 


[xni(k)-xB(k)}. 


(16) 


. .... 

.  Wksat-a: . .  iri -  ,  , 
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For  notational  convenience  we  define 

y,.  “  [xp(k)  +  v3(k  -  1  )T(k  *  1)  -  xm(k  +  1)]  (17) 

y  (18) 

Substituting  Eqs.  (15)  and  (16)  into  Eqs.  (13)  and  (14)  and  using  the  definitions  Eqs.  (17) 
and  (18),  we  find  that  Eqs.  (13)  and  (14)  are  identical  and  are  given  by 

r  py2y2T'k  +  in 

E[y1y2  +  <xy2y2+ — ^ - J=o.  (19) 

Hie  set  of  Eqs.  (13)  and  (14)  are  singular.  To  properly  constrain  the  gradients,  we  use  the 
relation  given  by  Eq.  (11)  as  /3  =  cr2/(2  -  a),  which  yields 


E 


c 


r 

a2 .  y2y2 


T(k  +  1)  1 


It*)  J 


-  y2y2  \  +  Of(2y2y2  -  y2y1  ]  +  2y1y2  =  0. 


(20) 


Finally,  since  the  statistics  of  y:  and  y2  are  not-  known,  the  expected  value  will  be  re¬ 
placed  with  a  time  average: 


“2  >2^2 


T(k  +  1) 

m 


~y 2^2  +  +  •=  o. 


(21) 


ADAPTIVE  FILTER  WITH  UNIFORM  UPDATES 

For  uniform  updates  T(k  +  1)  =  T{k)  =  1\  and  a  directly  relates  to  the  bandwidth  of 
the  filter.  Therefore,  Eq.  (21)  becomes 


-2y  i?2 

a  =  - 


(22) 


Defining  p, (ft)  =  yJ7£  and  p2(ft)  =  y2y2,  we  comPute  Pj(fc)  and  p2(ft)  by  passing  yjy2 
and  y2y2  through  low-pass  filters; 


px(fc)  =  3>i(*  - 1) +  (i  -^Kyay^) 


(23) 


p2(fe) s  3bp2(k  - 1)  (i  -  3"6)(y2y2)* 


(24) 


An  example  is  used  to  illustrate  the  filter’s  performance.  A  target  is  flown  away  from  the 
radar  at  a  speed  of  1000  ft/s,  it  makes  an  180-<legree,  3-g  turn,  and  then  flies  in  a 
straight  line  in  a  crossing  course  nenr  the  radar,  as  illustrated  in  Fig.  1.  The  radar  updates 
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Fig.  1  —  Target  trajectory  used  ir  examples 


the  track  every  4  s.  The  tracking  is  performed  in  range  and  the  range  measurement  error 
is  assumed  to  be  Gaussian  distributed  with  a  =  500  ft.  Constants  Ja  and  Jb  were  chosen 
to  be  0.819  and  0.91,  respectively.  The  filter’s  equations  are  summarized  below. 

Measure  xm  ( k  +  1) 

Pl(*)  "  W*  -  1)  +  (1  ~  4  ~  r>T  + 

P2(*>  =  - 1)  +  (1  -  y4.)l*m(*l  -*p(fe)l2 

f  -2?!^)  ] 

*  2p2(fe) -p1(fe)j 


fe  =  fe  +  l 


(25) 


\(kf 

p 

vt(k) 

L  T 


(1-0) 


"1)‘ 


»„•(*  -  1) 


l 

TJ 


t*m  (^)l  ■ 


Store  xp(fc),  vg(k  -  1),  and  xm  ih ) 
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xp(k  +  l)  =  xt(k)  +  Tvt(k) 

Repeat.  (25) 

The  track  vs  initialized  by  setting  all  positions  equal  to  the  first  measured  position  and  the 
velocity  equal  to  zero.  On  the  next  two  measured  positions  the  algorithm  is  operated 
normally,  but  we  constrain  a  and  0  to  be  equal  to  unity.  After  the  first  three  detections 
the  algorithm  is  operated  normally.  The  mean  values  of  the  bandwidth  a  and  the  error 
xq.(k)  -  xp(k)  are  plotted  as  a  function  of  time  in  Figs.  2  and  3.  We  find  that  the  band¬ 
width  decreases  when  there  is  only  straight-line  motion  c2  the  target  in  range  and  increases 
to  accommodate  the  target's  3-g  turn  and  accelerations  caused  by  the  target's  crossing 
course  in  front  of  the  radar.  The  mean  error  is  small  except  in  the  regions  of  high 
accelerations  '  i  which  the  filter  lags  the  target.  The  standard  deviation  of  the  error  de- 
creak's  as  the  filter’s  bandwidth  is  decreased  and  increases  when  the  filter’s  bandwidth  in- 
creaser,  which  occurs  to  accommodate  the  accelerations  of  the  target,  as  shown  in  Fig.  4. 
It  was  found  that  the  .  ystem  acted  much  better  when  the  gain  of  the  low-pass  filter  Eq. 
(23)  was  less  than  unity  (5T,  <  ). 

3-g  TURN 


TARGET  CROSSING  POINT 
I 


A 


TIME  (*) 

Fig.  2  —  Adaptive  adjustment  of  the  parameter  a  as  a  function  of  time 


The  adaptive  filter  was  compared  to  two  constant-bandwidth  filters  of  a  =  0.9  and 
0.2  respectively.  First,  looking  at  the  filter  of  a  =  0.2,  we  find  that  the  response  is 
sluggish  and  the  mean  errora  caused  by  the  high  accelerations  are  enormous,  'lowever, 
the  standard  deviation  does  settle  to  a  reasonably  small  value.  The  filter  with  a  =  0.9 
responds  rapidly  to  the  accelerations  because  of  its  wide  bandwidth,  but  the  standard 
deviation  of  the  error  does  not  achieve  a  small  value.  The  adaptive  system  adapts  to  the 
changing  environment  and  in  most  cases  yields  reasonably  small  errors. 


ADAPTIVE  FILTER  WITH  NONUNIFORM  UPDATES 

For  nonuniform  updates  a  does  not  directly  relate  to  the  natural  frequency  uQ,  but 
is  related  by 

(26) 


6 


STANDARD  DEVIATION  OF  XM(k)-Xf(k),  U 
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Fig.  3  —  Mean  error  as  a  function  of  time  for  constant  gain  and  adaptive  a  -  0  filter  operated 

with  constant  update  time 


Fig.  4  —  Standard  deviation  of  the  error  as  a  function  of  time  for  constant  gain  and  adaptive  filter 

operated  with  «'onstant  update  time 
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'Die  adjustment  of  the  bandwidth  co0  is  performed  by  solving  Eq.  (21)  using  the  defini¬ 
tion  Eq.  (26).  However,  the  time  averages  are  difficult  to  perform  and  yet  solve  for  oo0. 
We  introduce  an  approximation  for  Eq.  (26).  J 

2  to0T(k) 

“  1  +  2 tojQT(k)  *  <27) 

Substituting  Eq.  (27)  into  Eq.  (21),  one  finds 

A  — . — .  —  0.5C On 

iol[y2y2T<k)T(k  + 1)  +  y2y2T2(fe)  +  y^T2**)]  +  — - 


X&y2y2m+- 


“  2^  1  ■* 


0.5  , - , 

+  ~y  [vD =  °- 


The  time  averages  in  Eq.  (28)  are  computed  by  first-order,  low-pass  filters,  and  the  ad¬ 
justable  bandwidth  is  computed  by  solving  the  quadratic  equation.  The  target  trajectory 
described  in  Fig.  1  is  used  to  illustrate  the  filter  performance.  The  sampling  time  T(fe)  is 
assumed  to  be  a  uniformly  distributed  random  variable  that  can  take  on  values  of  0.5  to 
8  s.  The  filter’s  equations  are  summarized  below. 

g-  - 

=  e~03*nk*1) 

«i(fc> =  - 1)  +  <1  -  3&)[*p(fe)  +  vg(k  - 1  )T(k  +  l)~xm{k  + 1)1 

x  l*m(fc)-*p(*)] 

92<ft)  =  - 1) +  (1  -  ?6)[*p(fc)  +  vg(k  -  1  )T(k  +  1) -xm(k  +  1)] 

X  [*m(k)-xpmT{k) 

93(fe)  =  W*  -  1)  +  (1  -  3ft)[Xp(fc)  +  vt(k  - 1  )T(k  +  1) -xm(k  +  1)] 

X  ^m(k)-xp(k)]T(k)T(k) 

M*)  =  W*  -  1)  +  (1  ~  )'*„(*)  ~*p(*)l  l  *«(*)  -*p(*)JW 

?o (*)  =  - 1)  +  (1  -  5b)lxm{h)-xpm  [xm(k) -xp(k)]T(k)T(k) 

96(fe) '  ~  !)  +  tt  -  Jfr)[xm(fe)  -Xp(fe)]  [*m(fe)  -*p(*)J7(*)T(*  +  1) 


Wol7e<fc)  +  75<fe)  +  73(fe)l  +w0  f  $  [2q4(k)  +  3q2(k)J  + 


«»l-e 


-2£w0JT*n) 


8 
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P  (2 -a) 


k  =  k  +  l 


\l~a)  (1  - a)T(fe)~]  xg(k  -  1)1  ra"| 
-8  +  '  8 
T(k)  (1_/3  J  "*(fe "  1)J  Tjk) 


Store  xp(k),v,(k-l)  ,xm(k),7(k) 
xp(k  +  l)=xt(k)  +  T(k  +  l)vs(k) 

Repeat.  (29) 

Note  that  the  first  order  averaging  filter’s  time  constants  and  Jb  ai-e  adjusted  as  a 
function  of  the  sampling  interval  where  cjfl  =  0.05  and  u)b  =  0.024.  This  is  done  to 
maintain  a  reasonably  constant  bandwidth  in  time.  £  was  set  equal  to  0.75.  The  filter 
was  initiated  by  setting  all  positions  equal  to  the  first  measured  position  and  the  velocity 
equal  to  zero.  On  the  next  two  measurements  a  was  adjusted  by  setting  cuQ  =  0.2.  All 
subsequent  measurements  followed  the  algorithm.  Time  histories  of  the  bandwidth  and 
error  are  shown  in  Figs.  5  and  6  for  a  given  trial.  We  find  that  the  results  are  similar  to 
the  uniform  update  cases  (Figs.  2  and  3)  in  that  the  bandwidth  of  the  niter  is  wide  whe^, 
the  target  is  accelerating,  and  decreases  to  a  small  value  when  the  target  is  moving  in  a 
straight  line. 


N 

T> 

S  0.2 


3-g  TURN 


TARGET  CROSSING  POINT 


0  100  200  300  400 

TIME  (a) 

F  i$.  5  —  Adaptive  bandwidth  of  filter  aa  a  function  of  time  for  a 
random  update  time 
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Fig.  (J  —  Tims  history  of  error  for  adaptive  filter 
operated  with  a  random  update  time 


SUMMARY 

A  means  of  adaptively  adjusting  the  bandwidth  in  an  a  -  (3  filter  was  obtained  by 
using  the  LMS  criteria.  The  expected  value  operations  wore  approximated  with  time 
averages.  An  example  was  performed  to  show  that  the  bandwidth  of  the  filter  decreased 
when  the  target  was  flying  a  straight-line  course  and  opened  up  when  the  target  was  under 
acceleration  (real  and  apptrent).  The  filter  seemed  to  respond  reasonably  rapidly  to  its 
changing  environment  even  though  the  time  between  samples  was  quite  long. 
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